Hyperuniformity order metric of Barlow packings.

The concept of hyperuniformity has been a useful tool in the study of density fluctuations at large length scales in systems ranging across the natural and mathematical sciences. One can rank a large class of hyperuniform systems by their ability to suppress long-range density fluctuations through the use of a hyperuniformity order metric Λ[over ¯]. We apply this order metric to the Barlow packings, which are the infinitely degenerate densest packings of identical rigid spheres that are distinguished by their stacking geometries and include the commonly known fcc lattice and hcp crystal. The "stealthy stacking" theorem implies that these packings are all stealthy hyperuniform, a strong type of hyperuniformity, which involves the suppression of scattering up to a wave vector K. We describe the geometry of three classes of Barlow packings, two disordered classes and small-period packings. In addition, we compute a lower bound on K for all Barlow packings. We compute Λ[over ¯] for the aforementioned three classes of Barlow packings and find that, to a very good approximation, it is linear in the fraction of fcc-like clusters, taking values between those of least-ordered hcp and most-ordered fcc. This implies that the value of Λ[over ¯] of all Barlow packings is primarily controlled by the local cluster geometry. These results highlight the special nature of anisotropic stacking disorder, which provides impetus for future research on the development of anisotropic order metrics and hyperuniformity properties.

Inverse Design of Colloidal Crystals via Optimized Patchy Interactions.

Inverse statistical mechanics is a powerful optimization methodology that has been widely applied to design optimal isotropic pair interactions that robustly yield a broad spectrum of target many-particle configurations or physical properties. In this work, we generalize inverse techniques to design experimentally realizable spherical colloidal particles with optimized "patchy" anisotropic interactions for a wide class of targeted low-coordinated two-dimensional crystal structures that are defect-free. Our target crystals include square, honeycomb, kagomé, and parallelogrammic crystals. The square, honeycomb, and kagomé crystals possess desirable photonic, phononic, and magnetic properties, which are useful for a wide range of applications. We demonstrate that these target configurations can be robustly achieved with relatively few defects at sufficiently low temperatures. Our findings provide experimentalists with the optimal parameters to synthesize these crystals with patchy colloids under standard laboratory conditions.

Perspective: Basic understanding of condensed phases of matter via packing models.

Packing problems have been a source of fascination for millennia and their study has produced a rich literature that spans numerous disciplines. Investigations of hard-particle packing models have provided basic insights into the structure and bulk properties of condensed phases of matter, including low-temperature states (e.g., molecular and colloidal liquids, crystals, and glasses), multiphase heterogeneous media, granular media, and biological systems. The densest packings are of great interest in pure mathematics, including discrete geometry and number theory. This perspective reviews pertinent theoretical and computational literature concerning the equilibrium, metastable, and nonequilibrium packings of hard-particle packings in various Euclidean space dimensions. In the case of jammed packings, emphasis will be placed on the "geometric-structure" approach, which provides a powerful and unified means to quantitatively characterize individual packings via jamming categories and "order" maps. It incorporates extremal jammed states, including the densest packings, maximally random jammed states, and lowest-density jammed structures. Packings of identical spheres, spheres with a size distribution, and nonspherical particles are also surveyed. We close this review by identifying challenges and open questions for future research.

Can exotic disordered "stealthy" particle configurations tolerate arbitrarily large holes?

The probability of finding a spherical cavity or "hole" of arbitrarily large size in typical disordered many-particle systems in the infinite-system-size limit (e.g., equilibrium liquid states) is non-zero. Such "hole" statistics are intimately linked to the thermodynamic and nonequilibrium physical properties of the system. Disordered "stealthy" many-particle configurations in d-dimensional Euclidean space [Doublestruck R]d are exotic amorphous states of matter that lie between a liquid and crystal that prohibit single-scattering events for a range of wave vectors and possess no Bragg peaks [Torquato et al., Phys. Rev. X, 2015, 5, 021020]. In this paper, we provide strong numerical evidence that disordered stealthy configurations across the first three space dimensions cannot tolerate arbitrarily large holes in the infinite-system-size limit, i.e., the hole probability has compact support. This structural "rigidity" property apparently endows disordered stealthy systems with novel thermodynamic and physical properties, including desirable band-gap, optical and transport characteristics. We also determine the maximum hole size that any stealthy system can possess across the first three space dimensions.

Hyperuniformity variation with quasicrystal local isomorphism class.

Hyperuniformity is the suppression of long-wavelength density fluctuations, relative to typical structurally disordered systems. In this paper, we examine how the degree of hyperuniformity [[Formula: see text]] in quasicrystals depends on the local isomorphism class. By studying the continuum of pentagonal quasicrystal tilings obtained by direct projection from a five-dimensional hypercubic lattice, we find that [Formula: see text] is dominantly determined by the local distribution of vertex environments (e.g. as measured by Voronoi cells) but also exhibits a non-negligible dependence on the restorability. We show that the highest degree of hyperuniformity [smallest [Formula: see text]] corresponds to the Penrose local isomorphism class. The difference in the degree of hyperuniformity is expected to affect physical characteristics, such as transport properties.

Classical many-particle systems with unique disordered ground states.

Classical ground states (global energy-minimizing configurations) of many-particle systems are typically unique crystalline structures, implying zero enumeration entropy of distinct patterns (aside from trivial symmetry operations). By contrast, the few previously known disordered classical ground states of many-particle systems are all high-entropy (highly degenerate) states. Here we show computationally that our recently proposed "perfect-glass" many-particle model [Sci. Rep. 6, 36963 (2016)10.1038/srep36963] possesses disordered classical ground states with a zero entropy: a highly counterintuitive situation . For all of the system sizes, parameters, and space dimensions that we have numerically investigated, the disordered ground states are unique such that they can always be superposed onto each other or their mirror image. At low energies, the density of states obtained from simulations matches those calculated from the harmonic approximation near a single ground state, further confirming ground-state uniqueness. Our discovery provides singular examples in which entropy and disorder are at odds with one another. The zero-entropy ground states provide a unique perspective on the celebrated Kauzmann-entropy crisis in which the extrapolated entropy of a supercooled liquid drops below that of the crystal. We expect that our disordered unique patterns to be of value in fields beyond glass physics, including applications in cryptography as pseudorandom functions with tunable computational complexity.

The Perfect Glass Paradigm: Disordered Hyperuniform Glasses Down to Absolute Zero.

Rapid cooling of liquids below a certain temperature range can result in a transition to glassy states. The traditional understanding of glasses includes their thermodynamic metastability with respect to crystals. However, here we present specific examples of interactions that eliminate the possibilities of crystalline and quasicrystalline phases, while creating mechanically stable amorphous glasses down to absolute zero temperature. We show that this can be accomplished by introducing a new ideal state of matter called a "perfect glass". A perfect glass represents a soft-interaction analog of the maximally random jammed (MRJ) packings of hard particles. These latter states can be regarded as the epitome of a glass since they are out of equilibrium, maximally disordered, hyperuniform, mechanically rigid with infinite bulk and shear moduli, and can never crystallize due to configuration-space trapping. Our model perfect glass utilizes two-, three-, and four-body soft interactions while simultaneously retaining the salient attributes of the MRJ state. These models constitute a theoretical proof of concept for perfect glasses and broaden our fundamental understanding of glass physics. A novel feature of equilibrium systems of identical particles interacting with the perfect-glass potential at positive temperature is that they have a non-relativistic speed of sound that is infinite.

Transport, geometrical, and topological properties of stealthy disordered hyperuniform two-phase systems.

Disordered hyperuniform many-particle systems have attracted considerable recent attention, since they behave like crystals in the manner in which they suppress large-scale density fluctuations, and yet also resemble statistically isotropic liquids and glasses with no Bragg peaks. One important class of such systems is the classical ground states of "stealthy potentials." The degree of order of such ground states depends on a tuning parameter χ. Previous studies have shown that these ground-state point configurations can be counterintuitively disordered, infinitely degenerate, and endowed with novel physical properties (e.g., negative thermal expansion behavior). In this paper, we focus on the disordered regime (0 < χ < 1/2) in which there is no long-range order and control the degree of short-range order. We map these stealthy disordered hyperuniform point configurations to two-phase media by circumscribing each point with a possibly overlapping sphere of a common radius a: the "particle" and "void" phases are taken to be the space interior and exterior to the spheres, respectively. The hyperuniformity of such two-phase media depends on the sphere sizes: While it was previously analytically proven that the resulting two-phase media maintain hyperuniformity if spheres do not overlap, here we show numerically that they lose hyperuniformity whenever the spheres overlap. We study certain transport properties of these systems, including the effective diffusion coefficient of point particles diffusing in the void phase as well as static and time-dependent characteristics associated with diffusion-controlled reactions. Besides these effective transport properties, we also investigate several related structural properties, including pore-size functions, quantizer error, an order metric, and percolation thresholds. We show that these transport, geometrical, and topological properties of our two-phase media derived from decorated stealthy ground states are distinctly different from those of equilibrium hard-sphere systems and spatially uncorrelated overlapping spheres. As the extent of short-range order increases, stealthy disordered two-phase media can attain nearly maximal effective diffusion coefficients over a broad range of volume fractions while also maintaining isotropy, and therefore may have practical applications in situations where ease of transport is desirable. We also show that the percolation threshold and the order metric are positively correlated with each other, while both of them are negatively correlated with the quantizer error. In the highly disordered regime (χ → 0), stealthy point-particle configurations are weakly perturbed ideal gases. Nevertheless, reactants of diffusion-controlled reactions decay much faster in our two-phase media than in equilibrium hard-sphere systems of similar degrees of order, and hence indicate that the formation of large holes is strongly suppressed in the former systems.

Ground states of stealthy hyperuniform potentials: I. Entropically favored configurations.

Systems of particles interacting with "stealthy" pair potentials have been shown to possess infinitely degenerate disordered hyperuniform classical ground states with novel physical properties. Previous attempts to sample the infinitely degenerate ground states used energy minimization techniques, introducing algorithmic dependence that is artificial in nature. Recently, an ensemble theory of stealthy hyperuniform ground states was formulated to predict the structure and thermodynamics that was shown to be in excellent agreement with corresponding computer simulation results in the canonical ensemble (in the zero-temperature limit). In this paper, we provide details and justifications of the simulation procedure, which involves performing molecular dynamics simulations at sufficiently low temperatures and minimizing the energy of the snapshots for both the high-density disordered regime, where the theory applies, as well as lower densities. We also use numerical simulations to extend our study to the lower-density regime. We report results for the pair correlation functions, structure factors, and Voronoi cell statistics. In the high-density regime, we verify the theoretical ansatz that stealthy disordered ground states behave like "pseudo" disordered equilibrium hard-sphere systems in Fourier space. The pair statistics obey certain exact integral conditions with very high accuracy. These results show that as the density decreases from the high-density limit, the disordered ground states in the canonical ensemble are characterized by an increasing degree of short-range order and eventually the system undergoes a phase transition to crystalline ground states. In the crystalline regime (low densities), there exist aperiodic structures that are part of the ground-state manifold but yet are not entropically favored. We also provide numerical evidence suggesting that different forms of stealthy pair potentials produce the same ground-state ensemble in the zero-temperature limit. Our techniques may be applied to sample the zero-temperature limit of the canonical ensemble of other potentials with highly degenerate ground states.

Ground states of stealthy hyperuniform potentials. II. Stacked-slider phases.

Stealthy potentials, a family of long-range isotropic pair potentials, produce infinitely degenerate disordered ground states at high densities and crystalline ground states at low densities in d-dimensional Euclidean space R^{d}. In the previous paper in this series, we numerically studied the entropically favored ground states in the canonical ensemble in the zero-temperature limit across the first three Euclidean space dimensions. In this paper, we investigate using both numerical and theoretical techniques metastable stacked-slider phases, which are part of the ground-state manifold of stealthy potentials at densities in which crystal ground states are favored entropically. Our numerical results enable us to devise analytical models of this phase in two, three, and higher dimensions. Utilizing this model, we estimated the size of the feasible region in configuration space of the stacked-slider phase, finding it to be smaller than that of crystal structures in the infinite-system-size limit, which is consistent with our recent previous work. In two dimensions, we also determine exact expressions for the pair correlation function and structure factor of the analytical model of stacked-slider phases and analyze the connectedness of the ground-state manifold of stealthy potentials in this density regime. We demonstrate that stacked-slider phases are distinguishable states of matter; they are nonperiodic, statistically anisotropic structures that possess long-range orientational order but have zero shear modulus. We outline some possible future avenues of research to elucidate our understanding of this unusual phase of matter.

Confined disordered strictly jammed binary sphere packings.

Disordered jammed packings under confinement have received considerably less attention than their bulk counterparts and yet arise in a variety of practical situations. In this work, we study binary sphere packings that are confined between two parallel hard planes and generalize the Torquato-Jiao (TJ) sequential linear programming algorithm [Phys. Rev. E 82, 061302 (2010)] to obtain putative maximally random jammed (MRJ) packings that are exactly isostatic with high fidelity over a large range of plane separation distances H, small to large sphere radius ratio α, and small sphere relative concentration x. We find that packing characteristics can be substantially different from their bulk analogs, which is due to what we term "confinement frustration." Rattlers in confined packings are generally more prevalent than those in their bulk counterparts. We observe that packing fraction, rattler fraction, and degree of disorder of MRJ packings generally increase with H, though exceptions exist. Discontinuities in the packing characteristics as H varies in the vicinity of certain values of H are due to associated discontinuous transitions between different jammed states. When the plane separation distance is on the order of two large-sphere diameters or less, the packings exhibit salient two-dimensional features; when the plane separation distance exceeds about 30 large-sphere diameters, the packings approach three-dimensional bulk packings. As the size contrast increases (as α decreases), the rattler fraction dramatically increases due to what we call "size-disparity" frustration. We find that at intermediate α and when x is about 0.5 (50-50 mixture), the disorder of packings is maximized, as measured by an order metric ψ that is based on the number density fluctuations in the direction perpendicular to the hard walls. We also apply the local volume-fraction variance σ(τ)(2)(R) to characterize confined packings and find that these packings possess essentially the same level of hyperuniformity as their bulk counterparts. Our findings are generally relevant to confined packings that arise in biology (e.g., structural color in birds and insects) and may have implications for the creation of high-density powders and improved battery designs.

Precise algorithm to generate random sequential addition of hard hyperspheres at saturation.

The study of the packing of hard hyperspheres in d-dimensional Euclidean space R^{d} has been a topic of great interest in statistical mechanics and condensed matter theory. While the densest known packings are ordered in sufficiently low dimensions, it has been suggested that in sufficiently large dimensions, the densest packings might be disordered. The random sequential addition (RSA) time-dependent packing process, in which congruent hard hyperspheres are randomly and sequentially placed into a system without interparticle overlap, is a useful packing model to study disorder in high dimensions. Of particular interest is the infinite-time saturation limit in which the available space for another sphere tends to zero. However, the associated saturation density has been determined in all previous investigations by extrapolating the density results for nearly saturated configurations to the saturation limit, which necessarily introduces numerical uncertainties. We have refined an algorithm devised by us [S. Torquato, O. U. Uche, and F. H. Stillinger, Phys. Rev. E 74, 061308 (2006)] to generate RSA packings of identical hyperspheres. The improved algorithm produce such packings that are guaranteed to contain no available space in a large simulation box using finite computational time with heretofore unattained precision and across the widest range of dimensions (2≤d≤8). We have also calculated the packing and covering densities, pair correlation function g(2)(r), and structure factor S(k) of the saturated RSA configurations. As the space dimension increases, we find that pair correlations markedly diminish, consistent with a recently proposed "decorrelation" principle, and the degree of "hyperuniformity" (suppression of infinite-wavelength density fluctuations) increases. We have also calculated the void exclusion probability in order to compute the so-called quantizer error of the RSA packings, which is related to the second moment of inertia of the average Voronoi cell. Our algorithm is easily generalizable to generate saturated RSA packings of nonspherical particles.

Probing the limitations of isotropic pair potentials to produce ground-state structural extremes via inverse statistical mechanics.

Inverse statistical-mechanical methods have recently been employed to design optimized short-range radial (isotropic) pair potentials that robustly produce novel targeted classical ground-state many-particle configurations. The target structures considered in those studies were low-coordinated crystals with a high degree of symmetry. In this paper, we further test the fundamental limitations of radial pair potentials by targeting crystal structures with appreciably less symmetry, including those in which the particles have different local structural environments. These challenging target configurations demanded that we modify previous inverse optimization techniques. In particular, we first find local minima of a candidate enthalpy surface and determine the enthalpy difference ΔH between such inherent structures and the target structure. Then we determine the lowest positive eigenvalue λ(0) of the Hessian matrix of the enthalpy surface at the target configuration. Finally, we maximize λ(0)ΔH so that the target structure is both locally stable and globally stable with respect to the inherent structures. Using this modified optimization technique, we have designed short-range radial pair potentials that stabilize the two-dimensional kagome crystal, the rectangular kagome crystal, and rectangular lattices, as well as the three-dimensional structure of the CaF(2) crystal inhabited by a single-particle species. We verify our results by cooling liquid configurations to absolute zero temperature via simulated annealing and ensuring that such states have stable phonon spectra. Except for the rectangular kagome structure, all of the target structures can be stabilized with monotonic repulsive potentials. Our work demonstrates that single-component systems with short-range radial pair potentials can counterintuitively self-assemble into crystal ground states with low symmetry and different local structural environments. Finally, we present general principles that offer guidance in determining whether certain target structures can be achieved as ground states by radial pair potentials.

Effect of dimensionality on the percolation threshold of overlapping nonspherical hyperparticles.

We study the effect of dimensionality on the percolation threshold η(c) of identical overlapping nonspherical convex hyperparticles in d-dimensional Euclidean space R(d). This is done by formulating a scaling relation for η(c) that is based on a rigorous lower bound [Torquato, J. Chem. Phys. 136, 054106 (2012)] and a conjecture that hyperspheres provide the highest threshold, for any d, among all convex hyperparticle shapes (that are not a trivial affine transformation of a hypersphere). This scaling relation also exploits the recently discovered principle that low-dimensional continuum percolation behavior encodes high-dimensional information. We derive an explicit formula for the exclusion volume v(ex) of a hyperparticle of arbitrary shape in terms of its d-dimensional volume v, surface area s, and radius of mean curvature R[over ¯] (or, equivalently, mean width). These basic geometrical properties are computed for a wide variety of nonspherical hyperparticle shapes with random orientations across all dimensions, including, among other shapes, various polygons for d=2, Platonic solids, spherocylinders, parallepipeds, and zero-volume plates for d=3 and their appropriate generalizations for d≥4. Using this information, we compute the lower bound and scaling relation for η(c) for this comprehensive set of continuum percolation models across dimensions. We demonstrate that the scaling relation provides accurate upper-bound estimates of the threshold η(c) across dimensions and becomes increasingly accurate as the space dimension increases.

Microstructural degeneracy associated with a two-point correlation function and its information content.

A two-point correlation function provides a crucial yet an incomplete characterization of a microstructure because distinctly different microstructures may have the same correlation function. In an earlier Letter [Gommes, Jiao, and Torquato, Phys. Rev. Lett. 108, 080601 (2012)], we addressed the microstructural degeneracy question: What is the number of microstructures compatible with a specified correlation function? We computed this degeneracy, i.e., configurational entropy, in the framework of reconstruction methods, which enabled us to map the problem to the determination of ground-state degeneracies. Here, we provide a more comprehensive presentation of the methodology and analyses, as well as additional results. Since the configuration space of a reconstruction problem is a hypercube on which a Hamming distance is defined, we can calculate analytically the energy profile of any reconstruction problem, corresponding to the average energy of all microstructures at a given Hamming distance from a ground state. The steepness of the energy profile is a measure of the roughness of the energy landscape associated with the reconstruction problem, which can be used as a proxy for the ground-state degeneracy. The relationship between this roughness metric and the ground-state degeneracy is calibrated using a Monte Carlo algorithm for determining the ground-state degeneracy of a variety of microstructures, including realizations of hard disks and Poisson point processes at various densities as well as those with known degeneracies (e.g., single disks of various sizes and a particular crystalline microstructure). We show that our results can be expressed in terms of the information content of the two-point correlation functions. From this perspective, the a priori condition for a reconstruction to be accurate is that the information content, expressed in bits, should be comparable to the number of pixels in the unknown microstructure. We provide a formula to calculate the information content of any two-point correlation function, which makes our results broadly applicable to any field in which correlation functions are employed.

Effect of dimensionality on the continuum percolation of overlapping hyperspheres and hypercubes. II. Simulation results and analyses.

In the first paper of this series [S. Torquato, J. Chem. Phys. 136, 054106 (2012)], analytical results concerning the continuum percolation of overlapping hyperparticles in d-dimensional Euclidean space R(d) were obtained, including lower bounds on the percolation threshold. In the present investigation, we provide additional analytical results for certain cluster statistics, such as the concentration of k-mers and related quantities, and obtain an upper bound on the percolation threshold η(c). We utilize the tightest lower bound obtained in the first paper to formulate an efficient simulation method, called the rescaled-particle algorithm, to estimate continuum percolation properties across many space dimensions with heretofore unattained accuracy. This simulation procedure is applied to compute the threshold η(c) and associated mean number of overlaps per particle N(c) for both overlapping hyperspheres and oriented hypercubes for 3 ≤ d ≤ 11. These simulations results are compared to corresponding upper and lower bounds on these percolation properties. We find that the bounds converge to one another as the space dimension increases, but the lower bound provides an excellent estimate of η(c) and N(c), even for relatively low dimensions. We confirm a prediction of the first paper in this series that low-dimensional percolation properties encode high-dimensional information. We also show that the concentration of monomers dominate over concentration values for higher order clusters (dimers, trimers, etc.) as the space dimension becomes large. Finally, we provide accurate analytical estimates of the pair connectedness function and blocking function at their contact values for any d as a function of density.

Density of States for a specified correlation function and the energy landscape.

The degeneracy of two-phase disordered microstructures consistent with a specified correlation function is analyzed by mapping it to a ground-state degeneracy. We determine for the first time the associated density of states via a Monte Carlo algorithm. Our results are explained in terms of the roughness of an energy landscape, defined on a hypercubic configuration space. The use of a Hamming distance in this space enables us to define a roughness metric, which is calculated from the correlation function alone and related quantitatively to the structural degeneracy. This relation is validated for a wide variety of disordered structures.

Effect of dimensionality on the continuum percolation of overlapping hyperspheres and hypercubes.

We show analytically that the [0, 1], [1, 1], and [2, 1] Padé approximants of the mean cluster number S for both overlapping hyperspheres and overlapping oriented hypercubes are upper bounds on this quantity in any Euclidean dimension d. These results lead to lower bounds on the percolation threshold density η(c), which become progressively tighter as d increases and exact asymptotically as d → ∞, i.e., η(c) → 2(-d). Our analysis is aided by a certain remarkable duality between the equilibrium hard-hypersphere (hypercube) fluid system and the continuum percolation model of overlapping hyperspheres (hypercubes). Analogies between these two seemingly different problems are described. We also obtain Percus-Yevick-like approximations for the mean cluster number S in any dimension d that also become asymptotically exact as d → ∞. We infer that as the space dimension increases, finite-sized clusters become more ramified or "branch-like." These analytical estimates are used to assess simulation results for η(c) up to 20 dimensions in the case of hyperspheres and up to 15 dimensions in the case of hypercubes. Our analysis sheds light on the radius of convergence of the density expansion for S and naturally leads to an analytical approximation for η(c) that applies across all dimensions for both hyperspheres and oriented hypercubes. Finally, we describe the extension of our results to the case of overlapping particles of general anisotropic shape in d dimensions with a specified orientational probability distribution.